Light-harvesting in bacteria exploits a critical interplay 
between transport and trapping dynamics 
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Light-harvesting bacteria Rhodo spirillum Photometricum were recently found to adopt strikingly 
different architectures depending on illumination conditions. We present analytic and numerical 
calculations which explain this observation by quantifying a dynamical interplay between excita- 
tion transfer kinetics and reaction center cycling. High light-intensity membranes (HLIM) exploit 
dissipation as a photo-protective mechanism, thereby safeguarding a steady supply of chemical en- 
ergy, while low light-intensity membranes (LLIM) efficiently process unused illumination intensity 
by channelling it to open reaction centers. More generally, our analysis elucidates and quantifies 
the trade-offs in natural network design for solar energy conversion. 
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Photosynthesis is Nature's solution to the solar 
energy conversion problem [H-Q. Understanding what 
architectural designs it adopts, and why, arc important 
questions which could guide designs of future energy 
conversion devices. Many studies have clarified the 
exciton capture-transfer dynamics of reaction center 
(RC) pigment-protein complexes 0-0], their arrangement 
along the membranes that support them l3"15|, and 
even RC quantum effects [l6l. [l^ . Recent experimen- 
tal investigations resolved the locations of Light 
Harvesting (LH) complexes within the 2D membrane 
architecture of complete chromatophorc vesicles, to 
reveal an unexpected change in the ratio of complexes 
for bacteria grown under high (Fig. 1(a)) versus low 
illumination intensities (Fig. 1(b)). 

This Letter presents analytic and numerical results 
which explain this experimental observation [T^ by quan- 
tifying a trade-off which arises between two fundamental 
membrane requirements: (1) the need to convert large 
numbers of excitations into energetically useful charge 
separations within the RC, and hence promote metabolic 
activity, and (2) the need to avoid an oversupply of exci- 
tations and hence excessive bursts of energy, which could 
damage the photosynthetic machinery [lO, llH ■ Within 
our theory, the microscopic origin of this trade-off is the 
interplay between excitation transfer kinetics across the 
membrane architecture and reaction center cycling dy- 
namics. This generates a critical behavior of the mem- 
brane's efficiency when probed under different light in- 
tensities. Low light- intensity membranes (LLIMs) effi- 
ciently channel excess illumination intensity to open re- 
action centers, and hence are dominated by (1), while 
high light-intensity membranes (HLIMs) better exploit 
excitation loss through dissipation as a photo-protective 
mechanism in order to provide constant chemical energy. 



and hence are dominated by (2). Our analytic model pre- 
dicts a critical light intensity during growth, below which 
the synthesis of LH2 complexes should be dramatically 
enhanced. 

The photosynthesis process in purple bacteria 0, [Hj 
(see Fig. 1(c)) involves photons from sunlight being ab- 
sorbed with a rate ja = -\- ^2^2), where iVi(2) 
is the number of LHl (LH2) complexes, 71(2) 01 are the 
respective absorption rates per complex, and / is the 
light intensity. Each LHl contains one RC, i.e. the 
number of RCs is equal to A^i. Excitation transfer be- 
tween light harvesting complexes occurs at the follow- 
ing mean times (picoseconds) [1, 0, Il3|- ^12 = 15 for 
LH1-LH2; iai = 3.3 for LH2-LH1; tn = 20 for LHl- 
LHl; t22 = 10 for LH2-LH2; before reaching an RC at 
time ti^nc — 25 for LHl-RC. Once in an open RC, the 
special pair may become ionized {P^) on the timescale 
t+ = 3, eventually producing quinol (QBH2) by reducing 



quinone (Qs -^QbH) twice [ij. Otherwise back 

transfer occurs from RC to LHl with taci — 8. Be- 
fore a new Qb becomes available, the RC remains closed 
for a cycling time t during which an energetically useful 
charge separation is generated |19|. Dissipation 

through fluorescence or internal conversion happens at a 
constant rate 3 ■ The photo-excitation kinetics can 
be described by a collective population state vector p{t) 
which follows a master equation dtPi{t) = GijPj{t), 
where the element Gij of the rate matrix establishes the 
probability per unit time of a transition (due to absorp- 
tion, dissipation, RC ionization, transfer to neighbors) 
between collective states i and j. Due to small absorption 
rates, the probability that two excitations occupy either a 
single harvesting structure or RC is negligible. Light har- 
vesting complexes only have two possible states: no ex- 
citon present (unexcited) or one exciton state (excited). 
However the RC has four possible states: (un)excited 
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FIG. 1: (Color online) Top panels: Empirical architectures 
(i.e. digitized positions of LH complexes) from Ref. (a) 
High light-intensity membrane (HLIM), Jo = 100 W/m^ (b) 
Low light-intensity membrane (LLIM), Jo = 10 W/m^. Small 
orange dots: open reaction centers (RC) at snapshot during 
simulation. Large orange boxes: closed RCs. Large green 
circles: LHls. Small blue circles: LH2s. (c) Summary of 
dynamical processes of excitation transfer between LHl-LHl, 
LH2-LH2, LH1-LH2, LHl-RC. Dissipation (70), absorption 
(71(2)) and RC cycling time r (enlarged orange oval) are also 
shown. 



while being open or closed. Hence the state-space has 
size 2^i2^M^i = 2^^^+^\ 

Vesicles imaged with experimental Atomic Force Mi- 
croscopy (AFM) show that typically 2Ni A^2 ~ 
300 — 600. Given the large state-space, we study a 
discrete-time random walk simulation for excitations 
which includes dynamical coiipling to open/closed RC 
states, in contrast to Rcfs.jl, We use the empirical 
architectures to establish the likely neighbors (e.g. 
within 30A) between which excitations can hop, and im- 
plement the process of absorption and excitation transfer, 
dissipation, or RC capture (if within a RC) using Monte 
Carlo. All processes obey exponential distributions with 
mean values presented in Fig. 1(c). When two excita- 
tions reach a single RC, it closes for the cycling time r 



B Uj 3- We checked the accuracy of our stochas- 
tic numerical simulation by comparing its population- 
level predictions to those of a master equation for small 
chromatophores [l^l . The HLIM and LLIM architectures 
(Fig. 1, top panel) differ in the relative number of com- 
plexes and have stoichiometry s=N2/Ni. A typical snap- 
shot of open and closed RCs from our stochastic simula- 
tions (see Fig. 1 (a) and (b)) demonstrates that HLIM 
has fewer open RCs than LLIM in the experimen tally 



relevant regime of millisecond RC cycling times [18|, [19| . 
Our interest is the actual quinol output of the membrane, 
hence we calculate stationary-state obscrvables when the 
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FIG. 2: (Color online) (a) W{t) for LLIM (crosses) and 
HLIM (diamonds). The inset shows the quinol rate with re- 
spective RC cycling times for which LLIM has the same quinol 
output as HLIM, shown with arrows in the main plot, (b) Ef- 
ficiency 77 as a function of normalized intensity I /lo for t — 3 
ms in HLIM (diamonds) and LLIM (crosses) adapted mem- 
branes. Dotted lines are aids to the eye. In insets, distri- 
bution of open RCs, p{No), in HLIM (white bars) and LLIM 
(red bars) adapted membranes, are shown for the values high- 
lighted by arrows in main plots. Error bars smaller than the 
symbols are not shown. 



numerical simulations converge to a constant quinol rate. 
The quinol production rate W = ^ '^^Jt'^ ^^^^ ^^^^ 
at which nuc excitations produce ionization P'^. As- 
suming similar metabolic requirements under different 
illumination growth, the times suggested by arrows in 
Fig. 2(a) imply that LLIM has a shorter RC-cycling time 
than HLIM. This is consistent with greater quinone avail- 
ability in LHl clusters, as appropriate in LLIMs [ll| . 

Wc find qualitatively different behaviors in ly as a 
function of normalized intensity I / Iq (see Fig. 2(a) in- 
set). In HLIM, greater intensity does not change the 
quinol rate, while in LLIM higher illumination increases 
W. Therefore for higher light intensity, LLIM will be bet- 
ter than HLIM at processing potentially dangerous oc- 
currences of excess excitations. Due to fewer open RCs, 
HLIM will process only the necessary number of exci- 
tations for metabolism. The efficiency rj = nuc/uA is 
related to the quinol rate in the stationary limit through 
77 = 2W/jA, and quantifies the performance of a mem- 
brane in initiating RC ionizations from the ha total ab- 
sorbed excitations. Figure[2jb) shows that increased light 
intensity lowers 77 in both membranes due to a reduced 
number of open RCs (A^o); s^s shown by the distributions 
p{No) of open RCs at the top of Fig. 2(b). Consequently 
LLIMs have better efficiency than HLIMs since they have 
more open RCs in the high light intensity range, even 
though they have fewer RCs. 

An intriguing question arises as to whether clustering 
of LHls might help reduce the effective path that an exci- 
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differential equations: 



FIG. 3: (a) Ordered and (b) random membranes with s=8.09. 
(c) presents rj{T) for ordered (crosses) and random (diamonds) 
membranes. 



tation needs to take to a closed RC To explore this, 
we consider two extreme chromatophorc vesicles (Fig. 
|3lja) and (b)), both of which are compatible with relevant 
LLIM stoichiometries Their efficiencies are shown 

in Fig. ^c). In the experimentally relevant regime of 
millisecond RC cycling-time, open RCs are sparse - they 



are not clustered, and the architecture has no significant 
influence on 77. Less clustering does induce a slightly 
higher efficiency since the RC borders become easier to 
reach [20j . However apart from the benefit of quinone 
exclusion, clustering seems not to appreciably diminish 
the path length of excitations reaching a closed RC. On 
the other hand, clustering of LH2s might exclude active 
quinones from LHl-RC domains, thereby increasing their 
availability in core clusters and decreasing the RC-cycling 
time. Efficiency of the membranes will also depend on r, 
since this dictates how fast an RC opens due to a new 
quinone being available. If r is taken as a parameter, 
very small differences appear between membranes having 
equal stoichiometries but different network architectures. 
In Fig. 4(a), we compare the efficiencies of the represen- 
tative architectures presented in Fig. 1, calculated using 
the numerical stochastic simulations, as a function of RC- 
cycling time. In their respective illumination regimes, 
LLIM dissipates less excitations for r values in the bio- 
logically relevant millisecond range 3 l3| • This implies 
that even if their RC-cycling times are equal (i.e. no 
enhanced quinone diffusion in clustered LLIM), LLIM is 
efficient (77 ~ 85%) while HLIM provides a steady quinol 
supply by exploiting dissipation (77 w 20 — 40%). These 
findings hold irrespective of any other diffusion enhance- 
ments. 

Guided by these numerical findings, we now develop 
an analytic model to capture the underlying physics. As- 
sume that Ne excitations are created in the membrane 
at an absorption rate ja- Excitations either leave the 
membrane through dissipation at a rate , or at an RC 
at rate Ac where Ac depends on No- The number of 
RCs closing per unit of time is XcNe/2, while the num- 
ber opening per unit of time is -^(^^1 — No). Hence the 
number of absorbed photons is connected to the number 
of available (open) RCs by the following pair of coupled 



dt 
dNo 



-iXciNo)+7D)NE+"fA 



(1) 



hNi-No)-^Xc{No)NE (2) 

When Ac = 0, all RCs are closed. The maximum value 
(A^) occurs when all RCs are open. When the mem- 
brane is excited, the transfer-ionization rate per open 
RC is constant due to the fast excitation hopping rela- 
tive to the cycling time, i.e. Xc/Nq = X^/Ni which is 
also supported by numerical simulations (see Fig. 4(a) in- 
set). Given that 77 = X^NoNe /(Ni^a) , the steady state 
solution to Eqs. (|T]) and is pd^: 



V = 



{2Ni{X°c + id) + 1aX°ct - [m^iXl, + -fn? 



2iaXIt 

+4iVi7AA°,(7i, - XI)t + [lAXlrf] ^''~] 

In the limit of fast RC cycling-time (t— >0), 77 has the 
simple form 77 = (1 + 7£)/Ajl)^^. If all transfer paths are 
summarized by A^, this solution illustrates that 77 > 0.9 
, if the transfer-P reduction time is less than one tenth 
of the dissipation time in the absence of RC cycling. For 
finite r this analytic solution is in very good quantitative 
agreement with the numerical stochastic simulation, sup- 
porting our previously discussed interpretation (see Fig. 
mja)). Figure |4l[b) shows the complete analytical solution 
of Eqs. ([l]) and ([2]), in order to conffim the entire range 
of light intensities and RC closing times for which LLIM 
has a higher efficiency than HLIM. The assumed linear fit 
for Ac smears out an apparent power-law behavior. We 
have yet to find analytical solutions for 77 in cases where 
Xc{No) has a power- law dependence. 

Our analytical model can be used for easy compari- 
son of the metabolic outputs from experimentally dis- 
tinct AFM-imaged membranes, in order to provide addi- 
tional insight concerning the adopted and expected stoi- 
chiometries in Rsp. Photometricurr^^ . The bacteria are 
studied under different illumination conditions, assum- 
ing that comparable metabolic needs (i.e. quinol sup- 
ply) are accomplished in vesicles of area Aq . Our present 
aim is to find an expression for the quinol production 
rate W in terms of the environmental growth conditions 
and the responsiveness of purple bacteria through sto- 
ichiometry adaptation. In the stationary state, W = 
XcNe/'2 depends on the number of excitations within 
the membrane and on the details of transfer through 
the rate Ac- The LHl and LH2 complexes of area Ai 
and A2, respectively, fill a fraction p of the total vesi- 
cle area p = {AiNi + A2N2)/Ao. This surface occu- 
pancy has been shown to vary among adaptations, 
since LLIMs have a greater occupancy {p « 0.85) than 
HLIMs {p ~ 0.75) due to para-crystalline domains. The 
mean number of open RCs in the stationary state is 
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FIG. 4: (Color online) (a) Analytical model (continuous) and 
stochastic simulation (HLIM: diamonds, LLIM: crosses). In- 
set: Xc{No) from simulations within the linear portion of 
the main plot. Parameter values for LLIM and HLIM are 
A^={0.00771,0.0163}ps-^ and A^i = {40.01,70.66} respec- 
tively, consistent with empirical values, (b) ry as function of 
T and I/Io, obtained from complete analytical solution for 
LLIM (white) and HLIM (grey), (c) Numerical calculation 
of A5i,(s) (dots) vs. our analytic form from text (continuous), 
(d) W{s,I) as function of stoichiometry s and illumination 
intensity, with quinol rate contours of 1900 s~^ and 2100 s~^. 



^o = Ni- 2(j^lia) '^ = Ni- Wt. The linear Ac [No) as- 
sumption gives \c{s,W) = XUs) (1 - ^^I'tw^^^ ). 
The RC cycling time r(s) is expected to vary somewhat 
with adaptations due to quinone diffusion and different 
metabolic demands, and is described with a linear inter- 
polation using the values highlighted by arrows in Fig. 
^a). Likewise the rate Ac(s) must be zero when no RCs 
are present (s— >-oo), and takes a given value {to)~^ when 
the membrane comprises only LHls (s=0). Its depen- 
dence on s is satisfied by the form Xq{s) — (s/a+(io))~^, 
with adjustable parameter a, for several computer gener- 
ated membranes (see Fig. Sfc)). Solving Eqs.([T]) and ^ 
in the steady state, we obtain W — 2{\cis^wf+jD) '^^^'^^ 
can be solved to yield: 



2W{s,I) 




where B{s) = "^^^^^^^^^^^ . As can be seen in Fig. 4(d) in 
the high stoichiometry /high intensity regime, too many 
excitations would dangerously increase the cytoplasmic 



pH [H, H, 13 . Longer cycling times at higher light inten- 
sities are therefore helpful in order to keep power output 
bounded. The contours in Fig. Hl^d) of constant quinol 
production rate W, show that only in a very small in- 
tensity range will bacteria adopt stoichiometrics which 
are different from those experimentally observed in Rsp. 
Photometricum (sw4 and s«8) [loj . The empirical find- 
ing in Ref. [lol | that membranes with s=6 or s=2 are not 
observed, is consistent with our theory. More generally, 
our results predict a great sensitivity of stoichiometry 
ratios for 30-40 W/m^, below which membranes rapidly 
build up the number of antenna LH2 complexes. This 
prediction awaits future experimental verification. 

In summary, our analytic and numerical calculations 
elucidate and quantify the interplay which arises be- 
tween local (RC cycling) and extended dynamics (excita- 
tion transfer) in a chromatophore light-harvesting vesicle. 
In addition to explaining structural differences during 
growth, this new quantitative understanding may help 
accelerate development of novel solar micropanels mim- 
icking natural designs. 
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